Chemical Potential in the Gravity Dual of a 2+1 Dimensional System 
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We study probe D5 branes in D3 brane AdSs and AdSs-Schwarzschild backgrounds as a prototype 
dual description of strongly coupled 2+1 dimensional quasi-particles. We introduce a chemical 
potential through the \J(1)r symmetry group, U(l) baryon number and a U(l) of isospin in the 
multi-flavour case. We find the appropriate D5 embeddings in each case - the embeddings do not 
exhibit the spontaneous symmetry breaking that would be needed for a superconductor. The isospin 
chemical potential does induce the condensation of charged meson states. 



INTRODUCTION 



Recently there has been interest in whether the 
AdS/CFT Correspondence [J, S, Q can be used to un- 
derstand 2+1 dimensional condensed matter systems (for 
example 0, i, @, 0, S, @, E3] ) • The typical UV degrees of 
freedom in these systems are electrons in the presence of 
a Fermi surface and a gauged U(l), QED. When brought 
together in certain 2d states they can become relativistic 
and strongly coupled - possibly such systems might in- 
duce superconductivity too by breaking the gauge sym- 
metry. The philosophy, which may be overly naive, is to 
find relativistic strongly coupled systems that show these 
behaviours and hope they share some universality with 
the physical systems. Whether or not that linkage be- 
comes strong, it is interesting to study the AdS duals of 
2+ld systems. 

In this paper we will study the dynamics of the theory on 
the world volume of a mixed D3 and D5 brane construc- 
tion with a 2+1 dimensional intersection, which has pre- 
viously been studied at zero temperature in the absence 
of chemical potential in ll|, |l2|, |l3j . The gravity dual of 



the D3s, at zero temperature, is AdSs x S 5 , which is dual 
to the 3+1 dimensional Af — 4 super Yang Mills theory. 
Here these interactions will be used to loosely represent 
strongly coupled "phonons" . We will introduce 2+ld 
"quasi-particles" via D5 branes (with a 2+ld intersection 
with the D3s) - states connecting the two set of branes 
should be expected to carry quantum numbers that in- 
teract with the D3 brane dynamics and flavour quantum 
numbers associated with the number of D5 branes - the 
full field theory can be found in [HI]. We will work in 
the probe approximation for the D5 branes which corre- 
sponds to quenching quasi-particle loops in the phonon 
background At zero chemical potential the theory 

has Af = 4 supersymmetry and at zero quasi-particle 
mass is conformal [111 12, [l3|]. The system is related 
to the higher dimensional D3-D7 intersection where the 
Af = 4 gauge theory on the D3 branes has been used to 
describe gluon dynamics and the D3-D7 strings quarks - 
some progress in the study of the properties of mesons in 



3+ld strongly coupled gauge theories has been achieved 
[l5j ]. The D3-D5 defect system seems a natural starting 
point therefore for 2+1 dimensional systems. 

To attempt to mimic a solid state system one must 
weakly gauge a U(l) symmetry of the system and intro- 
duce an associated chemical potential (by setting A t = 
fi). There are a number of possible U(l)s that can play 
this role. 

Firstly the D3-D5 world volume theory has an unbroken 
SO (3) global symmetry corresponding to rotations in the 
4-plane transverse to the D5 brane. We will introduce a 
chemical potential for the quasi-particles with respect to 
a U(l) subgroup of the S0(3)- this can be done by simply 
spinning the D5 branes in an S0(2) plane [l6j]. The em- 
bedding of the D5 brane is described by a scalar that is 
charged under this U(l) symmetry so one naively expects 
to trigger superconductivity in the spirit described in [B[ 
- but here we would have an explicit understanding of 
the UV degrees of freedom the scalar describes. Naively 
one expects the scalar describing the D5 embedding to be 
destabilized by the presence of a chemical potential which 
gives the scalar a negative mass squared. An equivalent 
statement is that one expects the centrifugal force asso- 
ciated with the rotational motion of the brane to force 
it off the spin axis. In fact though we find the minimum 
area embedding for such spinning probe D5 branes and 
find this is not the case. 

The crucial physics is that the speed of light decreases as 
one moves into the centre of AdS - eventually it becomes 
less than the rotation speed of the D5 brane. We show, 
following the higher dimensional analysis in 1|| 17, 19| 
that there are regular D5 embeddings into the interior 
which have a more complicated embedding structure. 
The branes bend in the direction of the rotation so that 
there are two linked scalar fields describing the embed- 
ding - this richer theory turns out to not include super- 
conductivity, a subtlety on top of the arguments in [5[ . 

We can introduce mass terms for the quasi-particles that 
explicitly break the U(l) symmetry and we discuss the 
embeddings in these cases. There is a first order phase 
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transition when the R-charge chemical potential grows 
above the mass of the quasi-particle bound states - be- 
low the transition the quasi-particles exist as deconfincd 
particles whilst above it they are confined into bound 
states. This transition is analogous to the meson melting 
transition seen in this system and the D3-D7 system at 
finite temperature 2^, 21, 2^, 23[. We also analyze the 
finite temperature behaviour of these solutions by using 
the AdS§ Schwarzschild geometry as the background. 

Next we study a chemical potential for the U(l) associ- 
ated with baryon number for the quark fields - this seems 
the most natural candidate for how QED would manifest 
in the effective relativistic theory of a solid state system. 
The U(l) appears in the gravity dual as the U(l) gauge 
symmetry on the surface of the D5 branes - we allow con- 
figurations with non-zero profiles for these fields on the 
D5. Here we are again led by results in the D3-D7 system 
24|. 



At zero temperature the presence of the gauge field on 
the brane naively adds in an additional constraint on the 
solutions that that gauge field should be regular as the 
brane passes from positive to negative values of the effec- 
tive radial coordinate on the brane. This criterion rules 
out small perturbations of the standard flat embeddings 
of the D5 brane - instead the true solutions become ones 
where the D5 brane kinks through the origin of the space. 
The kink, which for large quark mass is rather sharp, has 
been interpreted in 24] as a tube of strings connecting 
the asymptotic branes. In j2(| it has been argued that an 
external charge could be responsible for the irregularity 
of the gauge field and that the flat embeddings should 
be retained. In either interpretation, amongst these con- 
figurations is one for zero quark mass which turns out 
to simply lie straight through the origin of the space - 
the chemical potential does not induce any R-charged 
operators to condense. The other fields on the D5 world 
volume carry no net baryon number and so do not couple 
to the chemical potential - there is no condensation. The 
system does not therefore act like a superconductor. 

The finite temperature behaviour of these solutions with 
non-zero baryon number chemical potential is also ex- 
plored. The D5 brane embeddings kink on to the black 
hole event horizon in this case. As in the D3-D7 case 
there is a first order phase transition between two differ- 
ent sorts of in-falling solutions which is analogous to the 
meson melting transition at finite temperature but zero 
chemical potential [H [HI, 0, . 



Finally we turn to isospin chemical potential in the case 
of multiple (but still probe) D5 branes. This theory seems 
less relevant to solid state physics because the quasi- 
particles have a U(2) flavour symmetry - electrons don't! 
On the other hand it is natural to discuss in this context 
and as advocated in may provide some lessons for p- 



wave superconductors. The embeddings of any individual 
brane is simply the same as for an equal magnitude bary- 
onic chemical potential and there is no induced R-charge 
breaking at zero quark mass. Where the theories differ is 
that there are vector bosons ('M /± ')on the branes' world 
volume that couple to the chemical potential - we work 
in a truncated version of the DBI action that is just Yang 
Mills theory on the D5 world volume. We show, with an 
analysis very similar to that of @, 0] and recent work in 
the D3-D7 system 27J, ]28| , that below some critical value 
of the temperature the W-bosons condense at a second 
order transition. This is dual to the formation of a spin 
one condensate which is charged under U(l) isospin num- 
ber in the gauge theory. Were one to identify that U(l) 
with QED we would have a superconductor. 



THE D3 THEORY 



We will represent the strong interaction dynamics with 
the large N M = 4 super Yang Mills theory on the surface 
of a stack of D3 branes. It is described at zero tempera- 
ture by AdS 5 x S 5 



ds 2 



TP- ax 3+l 



(1) 



+ lp^)(dp 2 + P 2 dn 2 2 + 



dr 2 



2 dn 2 ) 



where we have written the geometry to display the direc- 
tions the D3 lie in (£3+1 ), those we will embed the D5 
on (X2+1, p and ^2) and those transverse (r and f^). L 
is the AdS radius. 

At finite temperature the description is given by the AdS- 
Schwarzschild black hole 



2 r 2 

dg2 = ^{-h(u)dt 2 + dx 2 ) + -^tt-Ju 2 + L 2 dfl 2 (2) 



h{u) 



1 _ !fo 
u 4 



(3) 



It is helpful to make the change of variables to isotropic 
coordinates 



du dw 



Vu 4 - K w 



(4) 



and choose the integration constant such that if uq = 
the zero-temperature geometry is recovered 



2w 2 = u 2 + a/m 4 - ui 



(5) 
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The metric can now be written as 




+ £ (dp 2 + p 2 dn 2 2 + dr 2 + r 2 dttf) 

with w 2 = p 2 +r 2 , which shares the coordinate structure 
of ©. 



QUENCHED MATTER FROM A D5 PROBE AT 
T=0 



We will introduce quenched matter via a probe D5 brane. 
The underlying brane configuration is as follows: 

0123456789 

D3 

D5 • • • • 

In polar coordinates the D5 fills the radial direction of 
AdSs and is wrapped on a two sphere. 

The action for the D5 is just its world volume 



S ~ T J d 6 £V-detG ^ J dp pVl + r' 2 (7) 

where T is the tension and we have dropped angular fac- 
tors on the two-sphere. 

This is clearly minimized when r is constant so the D5 lies 
straight. The value of the constant is the size of the mass 
gap for the quasi-particles. We will mainly be interested 
in the conformal case where the constant is zero. Note 
the general large p solution is of the form 

r = m+- + .. (8) 
P 

Here m is an explicit mass term for the quasi-particles in 
the Lagrangian and c the expectation value for a bi-quasi- 
particle operator - note m has dimension one and c di- 
mension two adding to three as required for a Lagrangian 
term in 2+ld. The solution with non-zero c is not nor- 
malizable in pure AdSs. Note that when m = c = the 
theory is conformal. Including a non-zero more breaks 
the SO (3) symmetry ie it breaks one transverse SO (2) 
symmetry. From this it is apparent that m and c carry 
charge under that U(l). Were c to be non-zero when 
m = it would be an order parameter for the sponta- 
neous breaking of the U(l) symmetry. 



R-CHARGE CHEMICAL POTENTIAL/SPIN 

Our theory as yet lacks the relevant perturbation of the 
Fermi surface and the U(l) of QED. We will associate 
the U(l) with a subgroup of the SO(3) of the &2 - for 
concreteness we will use the angle in the X7 — xs direc- 
tions. 

To include a chemical potential we will spin the D5 brane 
in the angular direction cj> of this U(l) with angular speed 
p. 

The spinning of the D5 branes implies that the quasi- 
particles see a chemical potential. This is in fact a little 
bit of a peculiar limit since the background D3 theory also 
has fields, including scalars, charged under the U(l). We 
are not allowing that geometry to backreact to the chem- 
ical potential. In fact we had better not - the pure D3 
theory has a moduli space for separating the D3s in the 
transverse 6-plane. Were we to set them spinning they 
would scatter to infinity since there is no central force to 
support rotation. In the theory on the D3 surface there is 
a run away Bose-Einstein condensation. We simply wish 
to switch off this physics - it is not what we are inter- 
ested in - so we forbid such backreaction. The D3 theory 
is in an unstable state but will nevertheless provide some 
strongly coupled interactions for the quasi-particles that 
do see the chemical potential. 



An Overly Naive Ansatz 

We first look for solutions where the D5 embedding has 
4> = pt and we will allow the position r (the radial dis- 
tance in X7 — x%) to be a function of p. The action is 



S~ J dpp 2 ^{l + r' 2 )(l- {p2 L ^ r2)2 r 2 p 2 ) (9) 

Naively one is expecting the centrifugal force from the 
spinning to eject the brane from the axis at all but the 
end points where the boundary conditions hold the brane. 
This would lead to a spontaneous symmetry breaking or 
superconducting state. We will see that this is what this 
naive system tries to achieve. 

The equation of motion for r as a function of p is easily 
computed but unrevealing. At large p the solutions tend 

to the no-rotation limit r ~ m + - . 

p 

The (pair of) circle(s) in the (p, r) plane described by 
L 4 p 2 r 2 — (p 2 + r 2 ) 2 is clearly a zero of the action so 
branes wrapped there provide a solution to the equation 
of motion. Anything going within the locus described 
by the two circles is moving faster than the local speed 
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of light and is presumably not physical. This locus is a 
stationary limit surface - we call it the ergosurface below. 



There exist "Karch-Katz" type solutions [14j for D5- 
branes that do not encounter the stationary limit surface 

- these solutions essentially lie flat above everything plot- 
ted in Fig.l. We want to know what happens to those 
which have a close encounter with the ergosurface. It 
turns out that the curves which minimize the action like 
to hit the surface at a right angle. They then kink onto 
the surface where they can have zero action. 

The relation between m and c for curves impacting on 
the stationary limit surface in this way is shown in Fig.l 

- the presence of non-zero c at m = appears to indicate 
spontaneous breaking of the U(l) symmetry ie supercon- 
ductivity. Note there is a first order phase transition 
between the Karch-Katz embeddings and those hitting 
the ergosurface - we will discuss this transition further 
below. 




Ergosphere-impacting 



A More Sophisticated Ansatz 



We will now try a more sophisticated ansatz where the 
brane has in addition some profile 4>(w) where (f> is the 
angle on which they spin (ie <j> = /it + <fi(w)). The ansatz 



is inspired by the work in 17J where similar issues are 
encountered when a magnetic field is switched on on the 
brane's world-volume. 

We find it numerically convenient to switch coordinates 
and write the AdS geometry as 



J„2 _ w J 2 

as — jjax 3+1 



K(dw 2 

w z V 



d9 2 + sin 2 edtt 2 . + cos 2 Odtt 2 . 



(10) 



the D5 will now be embedded in the X2+\,w and £^2 
directions - the naive solutions above are recovered by 
looking for solutions that have 9(w) and <j> = fit where <f> 
is the 'first' angle of the Cl 2 - 

In these coordinates the Lagrangian for our more ambi- 
tious ansatz for the rotating D5 embedding is 



C = iu 2 sin 2 #x 



L-V 2 cos 2 



(1 + w 2 6' 2 ) + w 2 cos 2 6 <P' 2 



(11) 

If (j)' ~ fj, the two fi 2 terms compete against each other 
removing the naive intuition about centrifugal force. 

Since the action only depends on (f>' and not <f> one can 
integrate the equation of motion for cf>' . One could then 
substitute back in for </>' in terms of the integration con- 
stant - this though gives an action with a "zero over 
zero" form at the ergosurface that is hard to work with. 



Instead, 
£' = £- 



17, m 



following 

4>'§§7- This gives (setting 



we Legendre transform to 
§§r = J) 



FIG. 1: Embeddings of D5 branes impacting on the ergo- 
surface and some of the Karch-Katz type embeddings (note 
these actually exist down to the top of the ergosurface) (top). 
At the bottom is a plot of c vs m for embeddings impact- 
ing on the ergosurface (circles) and Karch-Katz embeddings 
(squares). The solutions oscillate around the value for the 
lowest Karch-Katz solution as the D5 approaches the very 
top of the ergosurface. 



The problem of course here is that the solutions are sin- 
gular at the ergosurface where they kink. This is a sign 
that our ansatz is wrong - none of this is the right physics. 



1 - 



L 4 /i 2 cos 2 ( 



„2a/2\ 



(12) 



6 sin 4 9 cos 2 0- J 2 ) 



This has a "zero times zero" form at the ergosurface 
which is simpler to work with numerically. 

For a solution that crosses the ergosurface we demand 
that the action be positive everywhere and this fixes J - 
the two terms must pass through zero and switch signs 
together. Having fixed J in this way one can then look at 
the 9 equation of motion near the ergosurface. Expanding 
near the surface, and after some algebra, one finds the 
following consistency equation for the 9 derivative 
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In the three-dimensional (w, 9, 4>) subspace the ergosur- 
face is the torus given by L 2 /icos9 = ±w, which in a 
plane of constant <f> gives two adjacent circles of radius 
^j-. Fig. 2 shows a sequence of regular solutions in the 
(p,r(p)) coordinates of the previous section. To obtain 
regular solutions one should make an odd continuation 
to the negative quadrant as shown. We show a full D5 
embedding in Fig. 3 with both the 9(w) and <fi(w) depen- 
dence plotted - note the D5 rotates at speed \i in the 4> 
direction (around the axis of the torus). 



FIG. 2: A selection of solution curves for D5 embeddings. 
The grey region is the interior of the ergosurface. 



Ergosphere- traversing 




FIG. 3: A particular solution curve in the three dimensional 
(w, 9, 4>) subspace. The torus represents the ergosurface. Note 
the D5 rotates at speed /i in the 4> direction (around the sym- 
metry axis). 



w 2 9' 2 + tan 9 w 9' - 1 = (13) 

There are thus two allowed gradients at the ergosurface. 
In fact numerically we find choosing any gradient focuses 
on to the same flow both within and outside the ergo- 
surface. We can numerically shoot in and out from a 
point near the ergosurface in order to generate regular 
embeddings. 



FIG. 4: A plot of c vs m for embeddings smoothly pene- 
trating the ergosurface (circles) and Karch-Katz embeddings 
(squares). The solutions again oscillate around the value for 
the lowest Karch-Katz solution as the D5 approaches the very 
top of the ergosurface indicating a first order transition. 



Clearly there is no spontaneous symmetry breaking in 
these solutions - the solutions smoothly map onto the 
solution which lies along the axis as the mass parame- 
ter m is taken to zero. In the field theory presumably 
the conformal symmetry breaking parameter (/x) which 
might trigger symmetry breaking is the same parameter 
as that telling us there's a plasma density cutting off the 
theory - there's no room for dynamics. This model turns 
out not to be an example of the behaviour studied in [f| . 

The presence of a non-trivial profile </>(w) for the em- 
beddings that penetrate the ergosurface indicates on the 
field theory side of the duality that there is a vev for the 
scalar field associated with the phase of the condensate c 
- this would be the Goldstone mode if there were sponta- 
neous symmetry breaking. Note that the regular Karch 
Katz embeddings, away from the ergosurface, have 4>(w) 
constant so there is no such vev. 

Again we see there is a first order transition between 
the Karch-Katz type solutions and those that enter the 
ergosurface region. We plot the values of c vs m for these 
solutions in Fig. 4 - it shows the same spiral structure 
around the first order transition as we saw with the naive 
ansatz. We will discuss the meaning of this transition 
below in the thermal context. 
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Thermal behaviour 



One can perform the same analysis in the thermal back- 

4 

ground. Writing 6 4 = there is again a torus- like 
ergosurface given by the equation 



L 2 \i cos 9 = ± 



1 w A -b A 



(14) 



and also a spherical horizon at w — b. One finds the 
horizon always lies within the ergosurface because the 
local speed of light is zero at the horizon. Note, below 
we find no phase transition when raising the temperature 
through the scale of the chemical potential. There would 
be a transition from a runaway Bose-Einstein conden- 
sation to a stable theory were we to allow the chemical 
potential to backreact on the geometry. 

One can form the Legendre-transformed Lagrangian 
(which recovers the T = case for b = 0) 



C 



1 w^+b 1 



WCQ W 



Y u> 4 (u> 4 +& 4 ) y 



L A y?(? B w 2 



jw' 



+b 4 ) 



^/(l + w 2 0' 2 )M 



„4„2 O 4 

b e c e 



-b 4 ) 2 



- J 2 - 



(15) 




FIG. 5: A selection of solution curves for D5 branes in the 
thermal geometry (with no = 1). The grey region is the 
interior of the ergosurface and the black region is the interior 
of the event horizon. 

The embeddings which extremize the action fall into two 
types - Karch-Katz type embeddings and those which 
hit the ergosurface. Fluctuations of the former would 
reveal a bound state spectrum. The latter embeddings 
inevitably fall onto the event horizon (a selection of these 
is plotted in Fig. 5 for uo = 1). In addition for these em- 
beddings that pass through the ergosurface gtt switches 
sign on the world volume - the ergosurface acts like a 



horizon for the world volume fields [191 ]. Here fluctua- 
tions would have a quasinormal spectrum along the lines 

of [13. 

There is therefore a first order transition in the behaviour 
of the theory as the quasi-particle mass goes through 
the scale of the chemical potential or temperature. This 
transformation is explored in detail in [l9[ . Note here it 
seems the transition is always a meson melting transition 
at finite temperature. At zero temperature the transi- 
tion is driven by quantum rather than thermal fluctua- 
tions and has been described in terms of a metal-insulator 
transition in 12511. 



The D3-D7 System 



Much of the above parallels results already found in the 
D3-D7 system 19]. That system describes an N = 2 
3+ld gauge theory with fundamental matter hypermul- 
tiplets in the gauge^ background of JV = 4 super- Yang 
Mills theory. In [161 ] an analysis similar to our "naive 
ansatz" was performed suggesting spontaneous symme- 
try breaking. Those authors have since refined their anal- 
ysis in a related system with a background electric field 
[1 71 ] and concluded that if regular embeddings are insisted 
upon the symmetry breaking is not present (see also [18j]). 
Were they to update [l6| they would find embeddings 
an alog ous to our D5 embeddings above as they indicate 
in [II]. 



BARYON NUMBER CHEMICAL POTENTIAL 



Another, and perhaps the most likely, way in which to 
embed the U(l) symmetry of QED into the brane set 
up is through the quasi-particle number global symme- 
try (essentially baryon number). The conserved vec- 
tor current and its source, which is effectively a back- 
ground gauge field configuration for this symmetry, man- 
ifest holographically as the U(l) gauge symmetry living 
on the world volume of the D5 brane. We can introduce 
a chemical potential for baryon number by switching on 
a constant A t component for this U(l) gauge field. We 
will study the embeddings of such a configuration at zero 
and non-zero temperature. Much of this analysis again 
mirrors that for the D3/D7 system which can be found 
in (24 1 . 
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Zero temperature 

The DBI action for the D5 brane including the surface 
gauge field is 



S ~ T 5 J <f^det(P[G ab ] + 2wa'F ab ) (16) 

We consider embeddings of the D5 brane in the p — r 
plane with in addition 2-na' Aq(p) — A(p) to represent 
the chemical potential. The action is then of the form 

C ~ p 2 y/l + r' 2 - A 12 (17) 




FIG. 7: Plot of parameters c, proportional to the condensate, 
versus m proportional to the quark mass in the case of baryon 
chemical potential for Q = 0.1. 



Since the action is independent of A the equation of mo- 
tion for A implies is a constant, Q. We find 



A' 2 = Q' 



1 + r' 2 
P A + Q 2 



(18) 



It is useful to perform a Legendre transform again (£' 
C - A'-^jt) and work with 



£' = v/(l + r' 2 )(p 4 + g 2 ) 



(19) 



The r independence again gives a simple equation for the 
embedding that is 



)2 _ r 2 



(20) 



with c\ a constant. 



FIG. 6: Embeddings with a baryon number chemical potential 
for Q = 0.1atT = with with nonzero ci. As ci approaches 
the numerical value of the charge Q the quark mass can be 
made arbitrarily large. 

To interpret this equation it is helpful to initially turn off 
the chemical potential, Q = 0. It is then clear that the 
solutions are singular at p — y/cy and the only regular 
case is c\ = so that r' — - we recover the usual flat 
embeddings. 

When we allow non-zero Q there become a bigger set of 
regular solutions - those with c\ < Q. These solutions 



are plotted in Fig. 6 for varying c\ and provide an alter- 
native embedding, that crosses through the origin, for 
each value of m in the large p asymptotics of the em- 
bedding. In 24j it was argued (in the D3-D7 case) that 
these are the true embeddings when there is a surface 
gauge field on the brane. We can see that the flat em- 
beddings (ci = 0) are not regular from (fl"8]) - they have a 
none zero gradient A' at p — so there will be a kink in 
the A field as it crosses over the r axis. For the solutions 
that pass through the origin though the A field is regu- 
lar. In [26] it has been argued that the irregular solutions 
should be maintained with the irregularity interpreted as 
the presence of an external source. 

As can be seen from Fig. 6, whichever interpretation is 
taken, the embedding for the case of massless quarks is 
unchanged from the usual flat embedding. The embed- 
ding does not therefore spontaneously break any symme- 
try with the introduction of a baryon chemical potential. 
Away from m = for the case of the regular embed- 
dings the embeddings do change and there is a conden- 
sate present. From ([20]) we can see that up to a sign ci is 
just the asymptotic parameter c that determines the con- 
densate. For small c\ the quark mass grows linearly but 
as Ci approaches Q m rises sharply - the resulting plot of 
c vs m therefore shows that the condensate asymptotes 
to a constant value for large mass - see Fig. 7. 

Another possible source of spontaneous breaking would 
be if the gauge field vev on the D5 led to other fields in 
the D5 brane world volume condensing. In fact though 
all the fields on the D5 are in the adjoint representation 
of, generically, a XJ(Nf) flavour symmetry. Adjoint fields 
of the U(l) of baryon number are chargeless and hence 
have no interaction with the gauge field. There is no 
possibility for such condensation and the system is not 
superconducting. 



8 



Finite temperature 

We can also study the theory with baryon number chem- 
ical potential at finite temperature by finding D5 embed- 
dings with a non-zero surface A t gauge field in the black 

4 

hole geometry (j6|). We again set b 4 = 

The DBI Lagrangian for such an embedding is 

a, A 



C 



+ b 4 / 1 {w 4 - 6 4 ) 2 



w 4 {w 4 + b 4 ) 

This time we have the gauge field 

l2 i f j£-W \ h , m 

(wi+b*) ) ( i+r ) 



(l + r> 2 )-A> 2 (21) 



A u = 



) +Q 2 



(22) 



The Legendre-transformed version of the Lagrangian is 



\ 



i (w 4 -b 4 y 

w 4 w 4 + b 4 



(1 + r' 2 ) p 4 



b 4 



(23) 



One can shoot out from the horizon attempting to fill out 
the m parameter space asymptotically. The solutions are 
very similar to those in the D3-D7 case as outlined in 
24j and see also 25| . There is no spontaneous symmetry 
breaking. There is a first order phase transition between 
the large m embeddings, that are essential flat except 
for a spike down onto the black hole, and smoother em- 
beddings that fall into the black hole at lower rn. This 
transition persists until the chemical potential becomes 
large, where there is no transition and the two phases 
coexist. In Fig. 9 we plot the quark condensate, c versus 
quark mass, m, for varying values of Q (which determines 
the chemical potential) around the critical value of Q 
where the phase transition between "spike" embeddings 
and smooth horizon entering embeddings ends. The dis- 
appearance of the phase transition is evident - the physics 
closely resembles that in the D3-D7 case discussed in de- 
tail in [24|. 



ISOSPIN CHEMICAL POTENTIAL 



The final possible source of a chemical potential in the 
D3-D5 set up is from the isospin symmetry present when 
there are two or more flavours of quasi-particle (D5) 
present. In contrast to the discussion of baryon num- 
ber above, there are clearly operators which carry isospin 
charge eg. {^oT^tj)) . These can be expected to condense 
at zero isospin chemical potential and break the symme- 
try spontaneously in the spirit of the zero temperature 




FIG. 8: Embeddings of the D5 branes for Q = 0.1 baryon 
chemical potential in units of the black hole temperature. 
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FIG. 9: Plot of the quark condensate, c versus quark mass, m, 
for varying values of Q (which determines the baryon number 
chemical potential) around the critical value of Q where the 
phase transition between "spike" embeddings (the top of the 
s-shape) and smooth horizon entering embeddings ends (lower 
part of the s-shape) . 



work in 29] and the phenomenological holographic model 
of p-wave condensation in @, Q . Here we will work first 
at finite temperature and only consider the case of zero 
quark mass - the embeddings of the branes are identical 
to those above for baryon number, where one uses the 
modulus of the isospin as the chemcial potential, so at 
zero quark mass the D5s lie straight along the axis as 
shown in Figure 6. The mesons of the theory are there- 
fore melted by the thermal bath but the operators can 
nevertheless condense. One is perhaps making a depar- 
ture from any obvious connection to a solid state system 
at this point since one would require a system with a 
U(2) or greater flavour symmetry on the quasi-particles 
- presumably there is only one sort of electron in a solid 
state system! 

The full theory of flavours on the D5 brane is expected 
to be unstable in the presence of a chemical potential at 
zero temperature. The situation is analogous to the D3- 
D7 system discussed in [301 ] - the squarks have a moduli 
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space in the gauge theory at zero temperature and chem- 
cial potential which shows up on the gravity side as a 
moduli space for the size of instanton configurations on 
the D7 (here D5s) world- volume [3JJ . An isospin chemi- 
cal potential will induce a negative mass squared for the 
scalars forcing the vev or instanton size to infinity. We 
will simply neglect this runaway behaviour here, fix the 
scalar vevs to zero and study the fermionic operators of 
the theory - hopefully this tells us about the behaviour 
of a theory with fermions but no scalars. 

The full DBI action to all orders in the surface gauge 
field is not fully known - an attempt to use the full DBI 
on a similar problem in the D3-D7 set up has recently 
appeared [27J . We though will follow the path of the D3- 
D7 analysis in [28| and just use the first order expansion 
of the action 



S ~ T 5 J d 6 £V-detG (\ - -Tr (F 2 )^j 



(24) 



We expect this action to represent the dynamics well. 

We will write the ansatz for the gauge fields as 

A = <P(p)r 3 dt + w{p)T 1 dx 1 (25) 

as in [8j . The coordinate representation of the Yang-Mills 
equation is 

dp (V^gF pin ) = -g Y M (S im 6 jl - 6 il S jm ) A^A^A vm 

(26) 



For our ansatz there are two equations of motion 

^ 5 °yv$ (27) 



-gg 00 g" '$')' = g YM , 



(\/=9<?V r O' = gYM^g^g 11 ® 2 ™ (28) 

Restricting ourselves to a massless quark D5 brane so 
the induced brane metric is given by isotropic AdS- 
Schwarzschild, the equations are (having absorbed fac- 
tors of gyM and L into the definitions of the fields) 



(r 4 + l)^ , 



f 



w $ 



— $ W 



(29) 
(30) 



In the isotropic coordinates one should shoot out from a 
small displacement x from the horizon using the initial 
condition w — wq and $ = $2^ 2 for constants uuq and 
$2. Near the boundary of AdS the solutions behave as 



<I> 



w 



P 

M 

r 

' £ 

r 



(31) 
(32) 



We search for solutions which have // = because these 
are solutions which are normalizable and hence describe 
a condensate of the charged mesonic operator. Requir- 
ing this to be the case one can solve the coupled non- 
linear equations to yield multiple branches of solutions. 
The lowest, monotonic branch is presumably the stable 
solution and for these solutions one can plot the depen- 
dence of the condensate on the chemical potential. Since 
the quarks we put in were massless, the only two scales 
are the chemical potential and the temperature and so 
large chemical potential can be equivalently viewed as 
low-temperature. 

The results are plotted in Fig. 10 - for \i <C T there is 
no condensation. For large fj, though there is a second 
order phase transition to a phase with a charged vector 
condensate. The behaviour is similar to that previously 
observed in @,@,|27|,|23|. As fi/T goes to infinity the con- 
densate c tends to a finite constant (« 0.3) when mea- 
sured in units of /x 2 , which is similar to the behaviour 
found in [28| for a pair of D7 probe branes - in their 3+1 
dimensional theory the condensate tends to a constant in 
units of fi 3 . 
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FIG. 10: Plot of the charged vector condensate (the parame- 
ter c from (J3TJ) ) in units of chemical potential squared versus 
in the case of isospin chemical potential at zero quark 
mass. 



SUMMARY 



We have proposed probe D5 branes in D3 brane back- 
grounds as a plausible dual for a strongly coupled quasi- 
particle theory in 2+1 dimensions - at zero tempera- 
ture and chemical potential the theory is supersymmet- 
ric and conformal. We introduced a chemical potential 
with respect to the global U(l) symmetries associated 
with R-charge and baryon number and found the result- 
ing regular D5 embeddings. These embeddings do not 
display spontaneous symmetry breaking and, indeed, at 
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zero temperature and zero intrinsic mass the theory is es- 
sentially indifferent to the chemical potential remaining 
as a state of conformal quasi-particles. 

For the R-charge case we show a first order phase tran- 
sition in the massive theory as the quasi-particle mass 
crosses the value of the chemical potential - on one side 
the quasi-particles are confined on the other they are not. 
At finite temperature the transition is between solutions 
that fall into the black hole and those that don't. 

At finite temperature in the baryon number case there is 
also a phase transition between two different black hole 
embeddings which is the equivalent of the meson melting 
transition at finite temperature but zero chemical poten- 
tial. If the chemical potential becomes too large then the 
transition ceases to occur for any quark mass. 

Finally we introduced isospin chemical potential in the 
case of two probe D5 branes and reproduced the second 
order transition to a phase with a charged vector conden- 
sate previously seen in other systems in 0, d, |27|, [H| 

We hope that these explorations will form a useful plat- 
form from which to find a holographic model of some real 
solid state system. We note that many of the transport 
properties of this system have also been recently explored 
in 13211 - 
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